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Abstract 

The spin-dependent structure functions g\(x), g%(x), gY W i x ) an d §2 (x) and 
their moments are studied in the CM bag model. The results show that 

(i) Jq 1 g2(x)dx = 0, i.e. the Burkhardt-Cottingham sum rule holds, hence 
g2 (x) must have at least one non-trivial zero besides x = and x = 1. 

(ii) Jq x 2 g2(x)dx is negative for the proton, neutron and deuteron. (iii) 
Jo x 2 g2(x)dx is about one order of magnitude smaller than Jq x 2 gi(x)dx, 
hence the twist-3 matrix element is approximately equal to the twist- 2 ma- 
trix element. The results are compared with most recent data and predictions 
from the MIT bag model, lattice QCD and QCD sum rules. 
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I. INTRODUCTION 



In recent years, the investigation of the spin structure of the nucleon by using deep 
inelastic lepton-hadron scattering (DIS) experiments has been an exciting and controver- 
sial field in hadron physics. The experiments performed by EMC, SMC at CERN []l]|§, 
and E142, E143 at SLAC |3|-|5| provide direct information on the matrix elements of spin- 
dependent operators in the nucleon. The spin-dependent DIS cross section is determined by 
the antisymmetric part of the hadronic tensor 

w£ = zw^VsiO^Q 2 ) + is* - p°^Mx,q 2 )} (i) 

where P and q are the four vectors of the nucleon and virtual photon momentums. 
gi(x,Q 2 ) and g 2 (x,Q 2 ) are two spin-dependent structure functions for the Bjorken vari- 
able x = Q 2 /2P ■ q = Q 2 /2Mu with Q 2 = —q 2 is the transfered four momentum squared 
and S a = U(P, S) r ) <J ^^U(P, S) is the covariant spin vector of the nucleon. 

According to Operator Product Expansion (OPE) analysis, to order M 2 /Q 2 , the lowest 
two moments of gi(x, Q 2 ) and g2{x,Q 2 ) can be written (QCD radiative corrections are not 
included) @-|nj 

£ 9l (x,Q 2 )dx = V)(Q 2 ) + ^ 2 \Q 2 ) + Ad^\Q 2 ) + Af 2 \Q 2 )] (2a) 

f 1 g 2 (x,Q 2 )dx = (26) 
Jo 

J Q x 2 9l (x,Q 2 )dx = - a ^(Q 2 ) + 0(—) (2c) 

j\ 2 g 2 (x,Q 2 )dx = -U 2 \Q 2 ) + \d^(Q 2 )+0(^) (2d) 

where a^°' 2 '(Q 2 ), d^°' 2 \Q 2 ) and f <y2 \Q 2 ) depend on the nucleon forward matrix elements of 
twist-2, twist-3 and twist-4 operators respectively. For example 

o(°>= £ e 2 Aq f , <P,S\^ flfl ^ f \P,S>=2A qf S, (3) 

qf=u,d,s... 

A-Qf {if = u,d,s, ..) are axial charges defined by the above axial- vectorial matrix elements. 
The singlet axial charge is proportional to the total quark helicity AS = Au + Ad + As in 
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the nucleoli. Using the semileptonic weak decay data, which are related to nonsinglet axial 
charges, the 1988 EMC result seemed to indicate that AE is surprisingly small and led to 
so-called "spin crisis" . Since then, an intensive study of g\ has been conducted. The efforts 
both from experimental and theoretical works led to a deeper understanding of internal spin 
structure of the nucleon, although many questions remain. The most recent reviews on this 
subject can be found in []Tl|— 113[| . Neglecting the terms of order M 2 /Q 2 in eqs.(2a-d), the 



longitudinal polarized structure function gi(x,Q 2 ) receives only twist-2 contributions. On 
the other hand, the structure function g 2 (x,Q 2 ) contains not only twist-2 but also twist- 
3 contributions corresponding to the matirx element d 2 (Q 2 ). The twist-3 contributions 
coming from spin-dependent quark gluon correlations do not vanish even in the large Q 2 
limit B,|!l-|16 . 



To separate the twist-2 and twist-3 contributions of g 2 (x, Q 2 ), one can write (Wandzura- 



Wilczek JT?j) 

g 2 (x,Q 2 )=g 2 vw (x,Q 2 )+g 2 (x,Q 2 ) (4) 
where the twist-2 piece (see, however, discussion in section III) is 

S 2 w (x, Q 2 ) = - 9l (x, Q 2 ) + f 1 %i(y, Q 2 ) (5a) 

Jx y 

and the twist-3 piece is 

Q 2 ) = 9i(x, Q 2 ) + g 2 (x, Q 2 ) - f 1 % 1 (y, Q 2 ) (56) 

Jx y 

It is easy to check that, to order 0(M 2 /Q 2 ), we have 

C 9Y W (x, Q 2 )dx = 0, f 1 x 2 gf w {x, Q 2 )dx = -\aW (6a) 
Jo Jo 6 

f 1 g 2 (x,Q 2 )dx = 0, f 1 x 2 g 2 (x,Q 2 )dx = ld (2) (66) 
Jo Jo o 

hence the twist-2 and twist-3 contributions of g 2 (x, Q 2 ) are separated. Measuring gi{x,Q 2 ) 
and g 2 (x,Q 2 ) as functions of x and Q 2 , one can obtain g 2 (x, Q 2 ), g 2 (x,Q 2 ) and their 
third and even higher moments. Up to order M 2 /Q 2 , g 2 (x, Q 2 ) uniquely measures the twist- 
3 contributions without involving any model-dependent analysis. Hence g 2 provides more 



detailed information about nucleon structure than does g\. Very recently, a preliminary 
experimental result of g 2 has been published |n| and more precise data obtained by E143 
experiment [pi at SLAC should be published soon. Since g 2 was discussed only briefly in a 



previous paper pt| , it is appropriate to give a detailed analysis in the modified Center-of- 
Mass (CM) bag model, using newly obtained data on g 2 . In section II, the CM bag model is 
briefly reviewed and the results for gz(x), g\ vw (x) and g%[x\ and their moments oS® ,2 \Q 2 ), 
S°' 2 \Q 2 ) are presented. Discussion and comparison with most recent data and other model 
results are given in section III. A brief summary is given in section IV. 

II. CM BAG MODEL 



As mentioned in a previous paper | 20|] , the bag model does not contain gluon fields 



explicitly, but the boundary of the bag-confined quarks simulates the binding effect coming 
from quark-gluon and gluon-gluon interactions (the gluon contribution to the proton spin 
in the bag model has been discussed by Jaffe [|ST| recently). Hence, the structure function 
g 2 calculated in the bag model does include higher twist effects. We have reported the CM 
bag model results |20| for the unpolarized and polarized structure functions F 1; F 2 , and gi, 
and briefly for g 2 . In this paper, we will focus our attention on g 2 . For reader's convenience, 
we give a brief review for CM bag model calculation of the structure functions. One needs 
to calculate the hadron tensor 

W, u (P,q,S) = ^J d i ye^<P,S\[J lx (y),J u (0)]\P,S> 

and separate the antisymmetric part to obtain Wfi v which can be expressed in general as 
Wt v = ie^vopicf I v)l a (x, Q 2 ) in the Bjorken limit, where I a (x, Q 2 ) depends on the model of 
the nucleon and the approximations used in the calculation. 

The basic assumptions and approximations of the CM bag model are: (i) the nucleon 
electromagnetic current J M (or 7NN vertex, see (2.7) in P0[]), can be approximately expressed 
by incoherent sum of single quark electromagnetic currents. It implies that the virtual 
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photon interacts with only one quark at a time and the other two quarks are spectators; 
this is an impulse approximation. The current includes not only the contribution of the 
struck quark but also those of the spectator quarks. Since the current satisfies translational 
invariance, four momentum is conserved, (ii) the nucleon consists of three valence quarks in 
their S'-wave state; higher excited states and higher Fock states which include gluons and 
sea quark pairs in addition to three valence quarks are neglected, (iii) SU(3) fiavor^S\J(2) spin 
wave functions for the proton and neutron are used. Symmetry-breaking effect is described 
in terms of a parameter £ = R% / R p d < 1, which simulates the smaller spatial size for the 



scalar u — d quark pair than that for the vector u — u and d — d quark pairs in the nucleon p2 
(iv) The effect of quark confinement due to nonperturbative quark-gluon and gluon-gluon 
interactions is described in terms of bound-state quark spatial wave functions, for instance 
the quark bag wave function in the cavity approximation in the MIT bag model or Gaussian- 
type quark wave function used in some other models. All necessary formulae for the CM bag 
model calculation can be found in [p0|] . A formal and general discussion on the theoretical 
basis of the CM bag model has been given in pi. We note that in addition to the CM bag 



model calculation, the transverse spin structure functions g 2 has also been computed in the 
original MIT bag model by Jaffe and Ji M, and other modified versions of the MIT bag 



model by Schreiber, Signal and Thomas (SST bag model [§4J]), and by Stratmann (MOD 
model H). 

Experimentally, gi(x,Q 2 ) and g 2 {x,Q 2 ) are measured by combining two different cases 
of deep inelastic scattering of polarized leptons on polarized nucleons: (i) the beam and 
target spin orientations are parallel, and (ii) the beam and target spin orientations are 
perpendicular. Experimentally, W^ v = ie^ ap (q p /v)I a (x,Q 2 ) can be calculated from various 
models of the nucleon. In this case it is convenient to choose suitable projection operators to 
extract g\ and g 2 from model results of I a (x, Q 2 ). One of possible projections is to extract 
gi = g L and g\+g 2 = 9t by choosing the nucleon spin parallel ('L') or perpendicular ('T') to 
the virtual photon momentum as we did in |]20|| . It should be noted that it is not necessary 
to choose the same projection as those used in the experimental analysis. 



Several parameters have been used in the CM bag model calculation. They are: (a) 'bag' 
radius R=5 GeV -1 , (b) SU(3) symmetry breaking parameter £ = 0.85 and (c) maximum 
momentum of quarks inside the nucleon |p max | = 0.6 GeV/c. R and £ were determined from 



the fit of the rms radius of the neutron and proton, and the ratio fM n / \i v |2!J . The model with 
these two parameters gives a fairly good result for the electromagnetic form factors of the 
nucleon and the magnetic moments of octet baryons. In particular, the neutron charge form 
factor is well reproduced. For the DIS parton distributions, the third parameter |p max | = 0.6 
GeV/c has been introduced. It constrains the unpolarized valence quark distributions to 
satisfy the sum rules 

/ u v (x)dx = 2, / d v (x)dx = 1 (7) 
Jo Jo 

As mentioned in [20|j, to compare the model results with data, the QCD evolution technique 
has to be used to evolve the parton distributions from the renormalization scale Ql to 
higher Q 2 where the experiments are performed. We choose <5o=0.81 (GeV/c) 2 and QCD 
scale parameter A=0.3 (GeV/c). The CM bag model with these parameters gives a good 
description for both unpolarized and polarized structure functions at 0.3 < x < 1, where 
the valence quark contributions dominate. For small- x region, < x < 0.3, the sea quark 
contributions are necessary. Using some QCD inspired phenomenological sea distributions 
(see (4.19) in p0[), we found that the first moment of gi(x,Q 2 ) is consistent with the 
experimental value and the first moment of gf(x, Q 2 ) — g"(x, Q 2 ) satisfies the Bjorken sum 
rule. For the higher moments (n > 2), the sea contributions coming mainly from the small x 
region are highly suppressed by the factor x n and are thus less important. Hence we neglect 
possible sea contributions for the higher moments of g 2 , gY an d gi- 

For QCD evolution of the structure function g\ and twist-2 piece g^ w from Q\ to 
Q 2 > Ql, the ordinary Altarelli-Parisi equations |26| can be used. For twist-3 part of 



g 2 , however, it has been shown [^7]] that due to mixing of twist-3 quark operators and 
quark-gluon operators with same twist and quantum numbers, the number of independent 
operators contributing to g 2 increases with n, where n refers to the n-th moment. It implies 
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that one cannot write down an Altarelli-Parisi type evolution equation for g 2 . This feature 



has been confirmed by several later calculations in p8|-|3~2|j and most recently in p3| . Hence 
there is no simple evolution equation for g 2 in the general case. One has to look for some 
approximate solutions. Two approximate evolution approaches under the limits N c — > oo 



or power n — > oo were suggested by Ali, Braun and Hiller P%fl . We use the approach in the 
large N c limit rather than the approximation in the large n limit, which only provides the 
asymptotic behavior of g2{x, Q 2 ) in the region x — > 1. We also note that as far as the third 
moments Jq 1 x 2 g 2 (x, Q 2 )dx are concerned, the Q 2 evolution is straightforward, i.e. a single 
power behavior of InQ 2 (for instance see ||33|| ) . 

In Fig. 1—5, we present the results of g{(x,Q 2 ), g 2 {x, Q 2 ), x 2 g 2 (x,Q 2 ), g 2 W (x, Q 2 ) 
(which is determined by gi(x,Q 2 )), and x 2 g 2 (x,Q 2 ) respectively. For the deuteron tar- 
get, gf(x,Q 2 ), x 2 g 2 (x,Q 2 ) and x 2 g 2 {x, Q 2 ) are shown in Fig. 6-8. All theoretical curves are 
calculated in the CM bag model at Q§=0.81 (GeV/c) 2 and evolved to Q 2 =5.0 (GeV/c) 2 
except for Fig. 4 and Fig. 6, where Q 2 =4.0 (GeV/c) 2 and Q 2 =3.0 (GeV/c) 2 respectively. 
Comparisons of our results for the third moments of gi and g 2 with recent data and other 
model predictions are listed in Tables I, II and III. The data for gi are taken from []T|,|2l,f3|,|3"4 
and those for g 2 are taken from [ITU 



III. RESULTS AND DISCUSSION 

1. For the leading term of first moments of gi(x), we obtain 

a proton = 0.252, ^neutron = —0.112, a deuteron = 0.064 (8) 

which can be compared with SMC data @ at < Q 2 >=5 (GeV/c) 2 : a { $ oton = 0.252 ±0.036, 
anlutron = -0.056 ± 0.024, a { ^ uteron = 0.046 ± 0.050 and E143 data |,§ at < Q 2 >=3 
(GeV/c) 2 : a$ oton = 0.254 ± 0.022, a^ eutron = -0.074 ± 0.027, af uterm = 0.084 ± 0.010. 

2. For the leading term of third moments of gi(x), we get 

^proton = 2-10 ■ 10" 2 , Onlutron = -1-86 • 10~ 3 , a deuteron = ^-74 ■ 10~ 3 (9) 
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while preliminary data at < Q 2 >=5 (GeV/c) 2 @ show: a$ oton = (2.42 ± 0.20) ■ 10~ 2 and 



''proton 

a< deuteron = (8-0 ± 1.6) • 10~ 3 . Comparisons with other models are listed in Tables I and III. 
We note that no g 2 data are available yet. 

3. Since g\ (x, Q 2 ) is always positive in the range < x < 1, hence Jq 1 x 2 g\[x)dx must be 
positive. From eq.(2c), the leading term of Jq 1 x 2 g\ (x)dx, i.e. a^J oton must also be positive 
as shown in eq.(9). However, the g±(x) as function of x is mostly negative except for large 

(2) 

x region, where g r {{x) is positive but very small. Hence its third moment, or a)J utron , is 
negative. For the deuteron, since the positive contribution from the proton is larger than 
the negative contribution from the neutron in the Jq 1 x 2 gf(x)dx, hence our model predicts a 
positive (ideuteron as shown in eq.(9) and Table III. 

4. As mentioned in section I, in the OPE approach including twist-2 and twist-3 operators 
in the presence of QCD corrections, the Burkhardt-Cottingham |3(| sum rule Jq 1 g 2 (x)dx = 
is known to hold. The CM bag model predicts 

f 1 g p 2 {x)dx = -0.0016, g%(x)dx = -0.0047 (10a) 

Jo Jo 

comparing to the numerical values for Jq 1 g[ p ' n ^dx in eq.(8), they are numerically consistent 
with zero. Hence in our model, the Burkhardt-Cottingham sum rule is satisfied within 
numerical errors. Most recent data |19| give 



f 1 gPJx) = -0.013 ± 0.028, f 1 g™(x)dx = -0.033 ± 0.082 (106) 

JO. 03 J 0.03 

which are consistent with zero. It should be noted that in the MIT bag model with SU(6) 
symmetry, g%(x) is identically zero by itself. However, the CM bag model with SU(6) 
symmetry breaking effects (£ ~ 0.85 < 1) predicts a nonzero g%(x). 

5. Since g 2 (x) is not identically zero in the model, the Burkhardt-Cottingham sum rule 
implies that g 2 (x) must change its sign at some x = xo, where < xq < 1, i.e. g 2 (x) must 
have at least one non-trivial zero. The question is whether the x-behavior of g 2 satisfies 

case 1 : g2(x) < 0, for x < Xq; ^(s) > 0, for x > Xq 



S 



or just opposite 

case 2 : 9i{x) > 0, for x < x ; ^(^) < 0, for x > x 

For case 1, one has Jq x 2 g2(x)dx > 0, while for case 2, Jq 1 x 2 g2(x)dx < 0. The CM bag model 
predictions for Jq 1 x 2 g2 > ' n,d \x, Q 2 )dx are all negative (see Table I, II and III or Fig.3 and 
Fig. 7). This implies the x-behavior of g2{x) fits case 2. The preliminary data [^| seem to 
favor our predictions (see, for instance, Fig.2 for the proton). 

6. For both the proton and deuteron, we now have a (2) > and d {2) - a (2) < 0. In the 
CM bag model, the moment Jq x 2 g2{x)dx (~ [d^ — a^]/3) is about one order of magnitude 
smaller than Jq 1 x 2 g\{x)dx (~ a^/2). It implies \d^ — a^\ << (recall \g 2 \ << gx), or 
^( 2 ) ~ a ( 2 )_ Hence the twist-3 matrix element d^ approximately equals to the twist-2 matrix 
element and the sign of d^ should also be positive for the proton and deuteron targets. 
This agrees with data |19[ but disagree with the negative sign predicted by the QCD sum 



rules p7| , |38[] and quenched lattice QCD pSfl . 

For the neutron, since Jq 1 x 2 g™(x)dx is negative and much larger than j Q l x 2 g2(x)dx in 

(2) 

magnitude, hence d nf ! utr . on is negative. This negative sign is consistent with the results given 
by other approaches P7|-|3T|. In magnitude, our result agrees with that given by the quenched 



lattice QCD, but much less than those given by the QCD sum rules. 

7. Our results show d^ ~ the twist-3 contribution is almost the same magnitude 

as twist-2 contribution. From (6a) and (6b), it implies that the g% (x) and #2 0*0 have 
opposite sign and approximately same magnitude. They almost cancel each other and 
lead to a very small g2{x) (see also Fig. 12a,b in ||20|| ). For the same reason the original 
Wandzura-Wilczek relation 



92(x,Q 2 ) = - gi (x,Q 2 )+ f 1 —gi(y,Q 2 ) 

J X v 



ix y 

is not a good approximation and the higher twist contributions may not be neglected. As 
pointed by Cortes, Pire and Ralston that the original Wandzura-Wilczek relation was 



derived by using the Dirac equation for free and massless quarks. Including quark mass 
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effect and gluon dependent term, an extended Wandzura-Wilczek relation was given (similar 
formula without quark mass term has been given by Jaffe J7|): 

g 2 ( X ,Q) = - gi ( X ,Q) + j x - gi (y,Q)--j x ~ Yy 

where the mass dependent term is another twist- 2 piece which is related to 'transversity' h?. 
The last term is the 'true' twist-3 piece arising from quark-gluon correlation. As emphasized 
in [4(| that when the quark mass term is included the twist-3 term cannot be isolated in a 
model independent way with a measurement of g± and g 2 . However, neglecting the strange 
quark contribution, the quark mass term (~ m q /M) should be negligible for up and down 
quarks. Hence the separation of twist- 2 and twist-3 pieces in eq. (4) seems to be a reasonable 
approximation. 

8. According to the OPE analysis, neglecting the M 2 /Q 2 corrections, the general for- 
mulae for the moments of the structure functions are 

f x n 9l (x, Q 2 ) = U n \Q 2 ), (n = 0, 2, 4, ...) (11a) 
Jo 2 



x n g 2 (x, Q 2 ) = -^_[aM(Q 2 ) - d (n) (Q 2 )] (n = 2, 4, 6, ...) (116) 



for n = 2, they reduced to (2c) and (2d). From (11a) and (lib), one obtains 

f x n [ 9l (x,Q 2 ) + ^±lg 2 ( x ,Q 2 )}dx = \d^ n \Q 2 ) (n = 2,4,6,...) (12) 
jo n 2 

hence one has for n = 2 

d< 2 ) = 2 f 1 x 2 [ 9l (x, Q 2 ) + \g 2 {x, Q 2 )]dx (13) 
Jo 2 

The CM bag model prediction for the function x 2 [gi(x, Q 2 ) + \g2{x,Q 2 )] for the proton is 
shown in Fig. 9 and that for the deuteron is shown in Fig. 10. One can see that in both 
cases, the model predictions for d*- 2 -* are nonzero and positive. This seems to be consistent 
with recent data fT9]j . 

9. If one assumes that (12) holds also for n = 1 one obtains 

[ 1 x[g 1 (x,Q 2 ) + 2g 2 (x,Q 2 )}dx = \d^{Q 2 ) (14a) 
Jo 2 
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To lowest order in a s , it was shown [|T^] by using the Field Theoretical Parton Model that 
c?W(Q 2 ) vanishes in the chiral limit, one has the Efremov-Leader-Teryaev (ELT) sum rule 

[ 1 x[g 1 (x,Q 2 ) + 2g 2 (x,Q 2 )]dx = (146) 

JO 

However, our results of g\ and g 2 do not satisfy this sum rule. To demonstrate this, we plot 
gf (x, Q 2 ) + 2g 2 J '' d \x, Q 2 ) as functions of x in Fig. 11 and Fig. 12 and compare them with 



recent proton and deuteron data |19| respectively. One can see that our model predictions 
are consistent with the SLAC E143 data. However, the ELT sum rule seems not to be 
supported by the data (at least for the proton data). 

IV. SUMMARY 

The study of transverse spin structure function g 2 (x, Q 2 ) has both theoretical and exper- 
imental interest. In the most naive parton model, the quark is asymptotic free and has no 
transverse momentum, g 2 (x) is identically zero. However, quarks inside the nucleon are not 
free, the binding effect, which arising from quark-gluon interactions, causes a nonzero trans- 
verse momentum for quarks and leads to g 2 ^ 0. Measurements of g 2 or g? = gi + g2 allow us 
to get more information about binding effects, which are mainly formulated as 'higher twist 
effects'. On the other hand, since Q 2 is large in the deep inelastic scattering, quark binding 
effects should not be significant and g 2 should be small, especially compared to g\. This 



seems to agree with most recent exprimental result Our model calculation is consistent 
with this conclusion. It should be noted, however, that the theoretical predictions given by 
different models or approaches seem not to fully agree with each other because of different 
approximations. In addition, as mentioned in section III that another twist-2 piece which 
is related to the quark mass and 'transversity' has been neglected in our discussion. On 
the experimental side, the errors of data are still quite large. We hope that several new 
experiments |4l^f43| to be performed in the next few years will provide more precise data 
and tell us more about the quark spin (including longitudinal and transverse) distributions 
in the nucleon. 
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Table I. Comparison of model results with proton data 
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Table II. Comparison of model results for the neutron 
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Table III. Comparison of model results with deuteron data 
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FIG. 2. CM bag model prediction for g%(x,Q 2 ) at Q 2 =5.0 (GeV/c) 2 , data from [19]. 
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FIG. 3. Same as Fig.2, but for x 2 g p 2 {x,Q 2 ) at Q 2 =5.0 (GeV/c) 2 , data from [19]. 
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FIG. 4. Evolution of twist-2 piece g% WW (x,Q 2 ) in the CM bag model at Q 2 =0.81 and 4.0 
(GeV/c) 2 , data from [19]. 
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FIG. 5. The twist-3 contribution x 2 g2(x,Q 2 ) calculated in the CM bag model and evolved to 
! =5.0 (GeV/c) 2 , data from [19]. 
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FIG. 6. The deuteron structure function gf(x, Q 2 ) calculated in the CM bag model and evolved 
to Q 2 =3.0 (GeV/c) 2 , data from [34]. 
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FIG. 7. The deuteron transverse structure function x^g^ix-, Q 2 ) calculated in the CM bag model 
and evolved to Q 2 =5.0 (GeV/c) 2 , data from [19]. 
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FIG. 9. CM bag model prediction for x 2 [gl{x,Q 2 ) + §#f(x,<2 2 )] at Q 2 =5.0 (GeV/c) 2 , data 
from [19]. 
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FIG. 10. Same as Fig.9 but for x 2 [gf{x,Q 2 ) + §c/f(x,Q 2 )] at Q 2 =5.0 (GeV/c) 2 . 
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FIG. 12. Same as Fig.ll but for the deuteron target, data from [19]. 
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